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We address the issue of bound state in the two-nucleon system in lattice QCD with the quenched 
approximation at the lattice spacing of a = 0.128 fm using a heavy quark mass corresponding to 
m K = 0.8 GeV. To distinguish a bound state from an attractive scattering state, we investigate the 
volume dependence of the energy difference between the ground state and the free two-nucleon 
state by changing the spatial extent of the lattice from 3. 1 fm to 12.3 fm. A finite energy difference 
left in the infinite spatial volume limit leads us to the conclusion that the measured ground states 
for not only spin triplet but also singlet channels are bounded. Furthermore the existence of the 
bound state is confirmed by investigating the properties of the energy for the first excited state 
obtained by 2x2 diagonalization method. The scattering lengths for both channels are evaluated 
by the finite volume formula derived by Liischer. 
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1. Introduction 

The strong interaction dynamically generates a hierarchical structure: three quarks are bound 
to form a nucleon with an energy of 1 GeV, and nucleons are in turn bound to form nuclei with 
a binding energy of 10 MeV or so per nucleon. This is a multi-scale physics that computational 
physics should explore, and lattice QCD is responsible for explaining the nature of nuclei based on 
first principles. 

Recently, we have made a first attempt to directly construct the Helium-3 and Helium-4 nuclei 
in quenched QCD [[I]] at a rather heavy quark mass corresponding to m n = 0.8 GeV, and success- 
fully confirmed the formation of Helium nuclei as a bound state. After our finding of the Helium 
nuclei, an evidence of the H di-baryon bound state in Nf ■ = 2 + 1 and Nf = 3 QCD were reported 
by NPLQCD [g, and HALQCD [§] Collaborations, respectively. The situation, however, is 
markedly different for deuteron, which is the simplest nucleus composed of two nucleons in the 
spin triplet channel, and yet evidence based on lattice QCD for bound state had never been reported 
before our paper [Q]. It is already quite some time ago that a first analysis of the two-nucleon system 
was made in quenched QCD & Much more recently, studies were made with a partially-quenched 
mixed action [Q] and Nf = 2 + 1 anisotropic Wilson action [§]. Extraction of the potential between 
two nucleons has been investigated in quenched QCD All these studies, however, tried to cal- 
culate the two-nucleon scattering lengths assuming, based primarily on model considerations with 
nuclear potentials, that the deuteron becomes unbound for the heavy quark mass, corresponding to 
m n ^ 0.3 GeV, employed in their simulations. 

To check the validity of this assumption, we need to investigate whether the bound state exists 
or not in the heavy quark mass region, where studies so far have been carried out, using the arsenal 
of methods solely within lattice QCD. If there is a bound state, the ground state energy never yield 
the scattering length if substituted into the Luscher's finite volume formula [[R]]. In such a case, 
the scattering length should be obtained from the energy of the first excited state. 

We carry out two types of calculations at a heavy quark mass corresponding to m n = 0.8 GeV 
in quenched QCD. The first one is a conventional analysis in which we investigate the volume 
dependence of the energy shift for the ground state. Different volume dependence is expected for 
scattering and bound states. In the second one we investigate the energy level of the first excited 



state employing the diagonalization method [11] to separate the first excited state from the ground 
state near the threshold of 2m^. If we find the ground state slightly below the threshold and the first 
excited state slightly above it, then such a configuration of the two lowest levels is consistent with 
the ground state being a bound state and the first excited state a scattering state with almost zero 
relative momentum. This method was previously used in a scalar QED simulation to distinguish a 
system with or without a bound state Jl2|]. 

Hereafter we call the analyses employed in the first and second calculations the single state 
and two state analyses, respectively. The results in this article have been reported in Ref. [|[]. 

2. Single state analysis 

We generate quenched configurations with the Iwasaki gauge action [ p4| ] at /3 = 2.416 whose 
lattice spacing is a = 0.128 fm, corresponding to aT x = 1.541 GeV, determined with ro = 0.49 fm 
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Figure 1: Effective energy shifts for 3 Si (left) and 'So (right) channels. Circle and square denote results for 
the source operators, respectively. 



as an input JT^]. We take three lattice sizes, L 3 x T = 24 3 x 64, 48 3 x 48 and 96 3 x 48, to investigate 
the spatial volume dependence of the energy difference between the two-nucleon ground state and 
twice the nucleon mass. The physical spatial extents are 3.1, 6.1 and 12.3 fm, respectively. 

We use the tadpole improved Wilson action with csw = 1-378 [15]. Since it becomes harder 
to obtain a reasonable signal-to-noise ratio at lighter quark masses for the multi-nucleon system, 
we employ a heavy quark mass at K = 0.13482 which gives m % = 0.8 GeV for the pion mass and 
niN = 1.6 GeV for the nucleon mass. Statistics is increased by repeating the measurement of the 
correlation functions with the source points in different time slices on each configuration. 

The quark propagators are solved with the periodic boundary condition in all the spatial and 
temporal directions using the exponentially smeared source q'(x,t) = I^Ae ~ B \* ~^q(y,t) after the 
Coulomb gauge fixing. On each volume we employ two sets of smearing parameters: (A,B) = 
(0.5,0.5), (0.5,0.1) for L = 24 and (0.5,0.5), (1.0,0.4) for L = 48 and 96. The onset of ground 
state can be confirmed by consistency of effective masses with different sources as shown later. 
Hereafter the nucleon operators using the first and the second smearing parameter sets are referred 
to as G\ and 62, respectively. 



2.1 Numerical results 

In order to determine the energy shift AEl = £aw — 2M^ precisely in each volume, we de- 
fine the ratio of the two-nucleon correlation function divided by the nucleon correlation function 
squared, R(t) = Gaw(?) / (Gjv(f)) 2 , where the same source operator is chosen for Gnn(0 and G^(f). 
The effective energy shift is extracted as A££ ff = In (R(t)/R (t + 1)). 

In the left panel of Fig. [I] we present typical results of time dependence of A£ , | ff for the 
sources in the 3 Si channel, both of which show negative values beyond the error bars in the plateau 
region of t = 8-11. Note that this plateau region is reasonably consistent with that for the effective 
mass of the two-nucleon correlation functions. The signals of AE[ S are lost beyond t 12 because 
of the large fluctuations in the two-nucleon correlation functions. We determine AEl by an expo- 
nential fit of the ratio in the plateau region, t = 8-12 for G\ and t = 7-12 for 62, respectively. We 
obtain a similar quality for the signal for the ! So channel on the (6.1 fm) 3 box shown in the right 
panel of Fig. [l[ 
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Figure 2: Volume dependence of energy shifts AEl as function of 1/L 3 for 3 Si (left) and 'So (right) chan- 
nels. The dashed line is extrapolation to the infinite volume, and filled circle denotes the extrapolated result. 



The volume dependence of the energy shift AEl for the 3 Si channel is plotted as a function 
of 1 /L 3 in the left panel of Fig. ^[ The results for the G\£ sources are consistent within the error 
bars. Little volume dependence for AEl indicates a bound state, rather than the 1/L 3 dependence 
expected for a scattering state, for the ground state in the 3 Si channel. The binding energy in 
the infinite spatial volume limit is extracted by a simultaneous fit of the data for the sources 
employing the fit function including a finite volume effect for the two-particle bound state [|2|, [TJ|], 




AE L = -—{ l + L ^ >, (2-1) 

m N 

where y and C r are free parameters, n is three-dimensional integer vector, and Yin denotes the 
summation without \n\ = 0. The binding energy, — AEoo, is determined from — AL^ = —y 2 /m^, 



assuming 2y mjj — J 2 — 2m^ m —y l jm^. The systematic error is estimated from the difference of 
the central values of the fit results choosing different fit ranges in the determination of AEl, and 
also using a constant fit as an alternative fit form. Adding the statistical and systematic errors by 
quadrature, we obtain — AE^ = 9.1(1.3) MeV for the binding energy. From the result, we conclude 
that the ground state in the 3 Si channel is a bound state. 

The right panel of Fig. || plots the volume dependence of the energy shift AL^ for the 
channel. Employing the same analysis as in the 3 Si channel, we find that — AEoo = 5.5(1.5) MeV 
in the infinite volume limit, which is 3.7 a away from zero. This tells us that the ground state in the 
l So channel is also bound at m n =0.8 MeV. Since the existence of the bound state in this channel 
is not expected at the physical quark mass, it might be a consequence of much heavier quark mass 
used in our calculation. 

3. Two-state analysis 

The focus of the analysis is the characteristic feature, well known from quantum mechanics, 
that the existence of a bound state implies a negative scattering length, and hence a scattering state 
just above the two-particle threshold in a finite volume. Our investigation is carried out with the 



diagonalization [11] of 2 x 2 correlation function matrix. 
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Figure 3: Effective energy shifts AE L a obtained from eigenvalues for ground (circle) and first excited 
(square) states in 3 Si (left) and (right) channels. Three solid lines and single line denote the result of the 
single state analysis with the G\ operator and free two-nucleon energy with the lowest momentum. 

We work with two spatial extents, 4.1 fm and 6.1 fm. The corresponding lattice sizes are 
L 3 x T = 32 3 x 48 and 48 3 x 48, respectively. The latter is the same size as in the first ensemble, 
but we regenerate independent configurations. Most of the simulation parameters, including the 
gauge and fermion actions, lattice spacing, quark mass, are identical to those explained in the 
previous section. The diagonalization method for the 2x2 matrix requires two operators each at 
source and sink time slice, which are explained in our full paper [Q]. We note that due to our criteria 
of the choice of the operators, this analysis does not provide an independent check for the ground 
state energy against the single state analysis. 

We diagonalize the following matrix at each t, M(t,to) = C(to)~ l C(t), where C(t) is the cor- 
relation function matrix of the two-nucleon operators and to a reference time. We determine the 
two eigenvalues X a (t) (a = 0, 1) of M(t,to) at each t and extract the energy of each eigenstate a 
through X a (t) = exp(—EL,a(t — to))- In order to determine the energy shift from the threshold as 
in the single state analysis, we define the ratio of the eigenvalue obtained from the diagonalization 
to the nucleon correlation function squared, R a (t) = A a (f)/(Gw(0) 2 - We also define the effective 
energy shift of the ratio R a as, AE^ a = \n(R a (t) /R a (t + 1)). 

3.1 Numerical results 

We show the results for the 3 Si channel on the (6.1 fm) 3 box. The effective energy shifts 
for the 3 Si channel are plotted in the left panel of Fig. ||. In Fig. || we use the nucleon correlator 
with the G\ operator for the denominator of R a (t). The ground state result AE^o * s reasonably 
consistent with the result of the single state analysis with the G\ source, which is expressed by the 
three solid lines in the figure. The first excited state is clearly higher than the ground state, but 
it is much lower than the free case with the lowest relative momentum, whose energy is given by 
2^J rnjj + (2n/L) 2 denoted by the single solid line in the figure. 

In the right panel of Fig. || the effective energy shifts for the ground and first excited states 
AE^lu for the 3 Si channel are shown. We find features similar to those in the 3 Si channel. We 
observe that the absolute value of the energy shift of the ground state is almost half of that in the 
3 Si channel. This is consistent with the observation in the first calculation. On the other hand, the 
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Figure 4: Volume dependence of energy shift AEl,i obtained from the eigenvalues of the first excited state 
for 3 Si and 'So channels as a function of 1 /L 3 . The statistical and systematic errors are added in quadrature. 



energy shift of the first excited state shown in the figure is almost twice larger than that in the 3 Si 
channel. This finding is consistent with the property of a system which contains a shallow bound 
state: The scattering length negatively increases as the binding energy decreases, diverging when 
the binding energy vanishes. From the results we confirm that the two-nucleon system in the both 
channels at the heavy quark mass of m n = 0.8 GeV has a bound state. 

In Fig. |] we plot the energy shift for the first excited state from the two lattice volumes as a 
function of 1/L 3 . A roughly linear behavior, with a larger shift on the (4.1 fm) 3 box compared 
to a smaller shift on the (6.1 fm) 3 box, is consistent with this state being a scattering state. We 
evaluate the scattering length using Luscher's finite volume formula [JTofa , where we find reasonable 
consistency between the two volumes [§]. If our finding of a bound state in quenched QCD at 
heavy quark mass smoothly continues to the physical point, then this is the first calculation which 
explained a negative scattering length for the deuteron channel. 



4. Conclusion and discussion 

We have carried out two calculations in quenched QCD to investigate whether the two nucleon 
systems have a bound state or not at the heavier quark mass, corresponding to m n = 0.8 GeV. In 
the first calculation, we have focused on the ground state of the two-nucleon system, and have 
investigated the volume dependence of the energy shifts obtained with two different source oper- 
ators. In the second calculation we have carried out two-state analysis using the diagonalization 
method with the 2x2 correlation function matrix. Based on these results we have concluded that 
the ground state is a bound state at the heavy quark mass in both the channels. 

While similar bound state in the 1 So channel was observed in recent Nf = 2 + \ QCD calcula- 
tion at m n = 390 MeV [§], the existence of the bound state looks odd from the experimental point 
of view. We expect that the bound state vanishes at some lighter quark mass, where the scattering 
length diverges changing the sign from negative to positive. Further reduction of the quark mass 
would decrease the scattering length. Confirmation of this scenario requires to investigate the quark 
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mass dependences of the binding energy and the scattering length. We leave this study to future 
work. 
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